Draft version November 21, 2009 

Preprint typeset using I^T^X style cmulateapj v. 08/22/09 



On 
O 

o 

> 

o 

(N 

o 

u 

43 

Oh- 
i 

o 

CO 

03 



> 

in 

ON 



On 
O 



13 



ARE NEWLY DISCOVERED HI HIGH VELOCITY CLOUDS MINIHALOS IN THE LOCAL GROUP? 

RlCCARDO GlOVANELLI 1 ' 2 , MARTHA P. HAYNES 1 ' 2 , BRIAN R. KENT 3 , ELIZABETH K. ADAMS 1 

Draft version November 21, 2009 

ABSTRACT 

A set of HI sources extracted from the north Galactic polar region by the ongoing ALFALFA 
survey has properties that are consistent with the interpretation that they are associated with isolated 
minihalos in the outskirts of the Local Group (LG). Unlike objects detected by previous surveys, such 
as the Compact High Velocity Clouds of Braun & Burton (1999), the HI clouds found by ALFALFA 
do not violate any structural requirements or halo scaling laws of the ACDM structure paradigm, nor 
would they have been detected by extant HI surveys of nearby galaxy groups other than the LG. At a 
distance of d Mpc, their HI masses range between 5 x 10 4 cP and 10 6 d 2 M© and their HI radii between 
< OAd and 1.6d kpc. If they are parts of gravitationally bound halos, the total masses would be on 
order of 10 8 -10 9 M©, their baryonic content would be signifcantly smaller than the cosmic fraction of 
0.16 and present in a ionized gas phase of mass well exceeding that of the neutral phase. This study 
does not however prove that the minihalo interpretation is unique. Among possible alternatives would 
be that the clouds are shreds of the Leading Arm of the Magellanic Stream. 

Subject headings: galaxies: spiral; — galaxies: distances and redshifts — galaxies: halos — galaxies: 
luminosity function, mass function — galaxies: photometry — radio lines: galaxies 



1. INTRODUCTION 

The ACDM paradigm which describes the evolution of 
structure predicts the existence of large numbers of low 
mass (< 10 9 Mq) halos. A cosmic census of dwarf galax- 
ies at z = indicates that such objects are rarer than 
expected from numerical simulations. In different guises, 
this circumstance has been referred to as the "the void 
phenomenon" (Peebles 2001) or the "missing satellite" 
problem (Klypin et al. 1999; Moore et al. 1999). The 
gap between the number of Milky Way (MW) satellites 
expected from numerical simulations and that obtained 
from observations has progressively narrowed in recent 
years with the discovery of fossil stellar structures in op- 
tical wide-field surveys. While this technique has de- 
tected numerous dwarf satellites of the MW, no gas-rich 
systems have been found and the number of known halos 
with circular velocities of v C i rc < 20 km s" 1 located in 
the general field remains very low. As is the case with 
the luminosity function of optical galaxies, the mass func- 
tion of extragalactic HI sources also has a faint end slope 
much shallower than predicted by the Press-Schechter 
formalism (Zwaan et al. 12006; Martin et al. 2009, in 
preparation). HI surveys are thus not detecting enough 
low mass sources to significantly alleviate the observed 
paucity of dwarf galaxies. 

Early on after their discovery, the possible extragalac- 
tic nature of HI High Velocity Clouds (HVCs) was con- 
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sidered (see previous work reviewed by Wakker & van 
Woerden 1997). More recently, the idea was forcefully 
revived by two groups: Blitz et al. (1999) and Braun & 
Burton (1999). The latter postulated that the so-called 
Compact HVCs (CHVCs), a subset of the HVC popula- 
tion defined by a relatively small angular extent and ap- 
parent spatial isolation are the baryonic counterparts of 
isolated, low mass halos distributed across the LG. Here- 
after we shall refer to this idea as the "minihalo/HVC 
hypothesis" . Two important problems mire this idea. 
The first is that, if present in other nearby groups of 
galaxies, CHVCs could have been detected by extant HI 
surveys; they have not (e.g. Pisano et al. 2007). The 
second is that the structural properties of the CHVCs 
are inconsistent with those expected for low mass halos, 
according to the ACDM structure formation scenario, as 
pointed out by Sternberg, McKee & Wolfire (2002, here- 
after SMW02): as minihalos, CHVCs appear undercon- 
centrated and, at typical LG distances — on order of 1 
Mpc — they would be too large. 

A widely adopted explanation of the discrepancy be- 
tween the numbers of observed and predicted dwarf 
galaxies is that gas accretion onto low mass halos is sup- 
pressed by the intergalactic UV radiation field (Ikeuchi 
1986; Rees 1986; cf. Hoeft et al. 2006). In their Mare 
Nostrum simulations, the latter show that the baryon to 
total mass fraction may fall steeply from the cosmic value 
of ~ 0.16 to only a few percent between halo masses 10 10 
M Q and M < 10 9 M . SMW02 showed that most of the 
gas in low mass halos should be found in a thermally 
stable, ionized phase which envelops a warm but neutral 
and stable component of much lower mass. Within the 
latter, a cold core may be able to form and partially con- 
vert into stars. A 10 9 M© halo may thus contain fewer 
than 10 6 Mq in a phase detectable optically or via its 
21 cm emission, thus explaining the paucity of such ob- 
jects in wide field surveys. The simulations by Ricotti 
(2009) improve somewhat the observational prospects: 
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he points out that, as the Universe expands, the inter- 
galactic medium cools and the UV radiaton background 
is diluted; changes in the IGM Jeans mass and the halo 
mass concentration then may make it possible for a previ- 
ously 'dormant' low mass halo to resume accreting gas at 
low z. He proposes that the recently discovered, nearby 
dwarf galaxy Leo T (Irwin et al. 2007) may be such an 
object. At a distance of 420 kpc, Leo T is a star forming 
galaxy with an HI mass of 2.8 x 10 5 M Q , an HI radius 
of 300 pc, an indicative dynamical mass within the HI 
radius of <~ 3.3 x 10 6 M Q , a total mass to V-band lu- 
minosity within the HI radius of 56 and a stellar mass 
of - 1.2 x 10 5 M Q (Ryan-Weber et al. 2008; see also 
Grcevich & Putman 2009). 

The precipitous drop in the cold baryon fraction for ha- 
los of mass < 10 10 M Q proposed by the simulations by 
Hoeft et al. and the models of SMW02 can explain the 
shallow slopes of the low end of optical and HI luminos- 
ity functions. It also indicates that the minihalo/HVC 
hypothesis could only apply to sources of more extreme 
properties than previously thought: HI masses of < 10 6 
Mq, sizes of < 1 kpc and linewidths of < 30 km s" 1 . 
The requirements for detection of such sources are very 
challenging, exceeding the combination of sensitivity and 
spectral resolution of most large-scale extragalactic HI 
surveys. However, the currently ongoing Arecibo Legacy 
extragalactic HI survey (ALFALFA) offers an opportu- 
nity for detection of systems matching the requirements 
mentioned above, albeit only to distances of a few Mpc. 
The question motivating this paper is whether candidates 
for minihalos with a detectable HI component exist. If 
they do, the safest place to look for them is away from 
a massive galaxy like the MW. We found them, but we 
still treat the minihalo interpretation with skepticism un- 
til the detection is achieved of (i) a stellar counterpart 
or (ii) of an internal velocity field that reveals more reli- 
ably a dynamical mass than a plain lincwidth or (iii) of 
analogous objects in other galaxy groups. 

2. OBSERVATIONS 

We report the results of 21cm HI observations, made 
with the Arecibo 305m radio telescope as part of the AL- 
FALFA survey. They cover the region between 9.5 h < 
RA < I6.5 h , +4° < Dec < +16°, of -1300 square de- 
grees, with a spectral resolution of — 5.5 km s _1 , an 
angular resolution of ~ 3'. 5 and a sensitivity such that 
an HI mass of ~ 5 x 10 4 M© can be detected at 1 Mpc 
distance. A full description of the observational mode of 
ALFALFA is given in Giovanelli et al. (2007), while the 
definition and goals of the survey are described in Gio- 
vanelli et al. (2005). Aspects of the processing pipeline 
are described in Saintonge (2007) and Kent (2008). 

The chosen region of the sky is a strip extending be- 
tween b — 35° and b — 80° in the North Galactic po- 
lar cap, well suited to minimizing the kinematical intru- 
sion of Galactic HI. Within that region, we searched for 
sources not associated with known galaxies with the fol- 
lowing additional criteria: (i) heliocentric velocity V© > 
120 km s _1 or Vq < —110 km s _1 (the V© range between 
-110 and +120 km s _1 was excluded in order to avoid 
confusion with Galactic emission); (ii) major HI diameter 
< 15'; (iii) no visible low surface brightness connection to 
features in the -110 to +120 km s _1 range. Twentyone 
ultra-compact clouds (hereafter UCHVCs) with positive 



TABLE 1 
Cloud Properties 



Property 



Avg. 



Range 



Vq, positive, km s 1 




124:320 


Vq, negative, km s _1 




-118:-142 


W 


24 km s" 1 


16:56 


R HI d lLc 


0.75 kpc 


,0.4:1.6 


5.25 


4.75:5.97 


log N HI [crrr 2 } 


19.10 


18.80:19.58 




7.35 


6.83:7:77 


t d 1 

Tcrossa Mpc 


0.08 Gyr 


0.02:0.32 



velocities (6 are unresolved) and 6 clouds with negative 
velocities were found. All sources detected by ALFALFA 
with a signal-to-noise S/N < 10 were re — observed and 
confirmed. Figure 1 shows an integrated flux map and 
a line profile of one of the detections. A catalog with a 
detailed description of individual sources is presented in 
Giovanelli et al. (2009, in preparation). Table 1 displays 
the mean and range values for the properties of the 27 
clouds, scaled by the unknown distance in Mpc, namely: 
the range in heliocentric velocities for positive and nega- 
tive velocity clouds, in km s _1 ; the full linewidth at half 
power W, of generally Gaussian shape; the mean radius 
at the isophote encircling 50% of the flux Rhi in kpc; 
the log 10 of the HI mass and of the total mass within 
Rhi, in solar units; the crossing time T cross = 2Rhi/W 
in Gyr. On the assumption that the clouds are parts of 
self-gravitating systems, the total mass is estimated via 



M tot (< R H i) ~ R H i<J 2 /G 



(1) 



with a — VF/2\/21n2. The "minimum intrusion" ap- 
proach (Giovanelli et al. 2005) adopted for data taking 
by ALFALFA delivers exceptional bandpass stability and 
allows reliably recovery of flux and size estimation for the 
UCHVCs. The possibility of extended envelopes of dif- 
fuse gas below the column density limit of the survey 
cannot however be excluded. 

The sky distribution of the UCHVCs is shown in Fig. 
2. Also plotted are the locations of galaxies with pri- 
mary distances of < 2.6 Mpc, as listed in the Catalog of 
Nearby Galaxies of Karachentsev et al. (2004). Those 
galaxies are, respectively from West to East: D634-03, 
Leo T, Sex B, Leo I, GR8, KKH86 and DDO 187. All 
but Leo I (the nearest at d = 0.25 Mpc) are detected 
in HI. Leo T and DDO 187 are just outside the Dec. 
range of the HI clouds. Figure 3 shows the location of 
clouds and nearby galaxies in the Galactic longitude vs. 
velocity in the Galactic Standard of Rest plane, where 
Vgsr = Vi sr + 225 sin I cos b and for Vi sr we assume a so- 
lar motion of 20 km s" 1 towards I = 57°, b = 25°. The 
three negative velocity clouds near RA = 16 h could pos- 
sibly be associated with the extended, perigalactic HVC 
complexes A and M. 

3. DISCUSSION 

The sky and velocity distributions of the clouds found 
by ALFALFA matches well that of dwarf galaxies in the 
LG. If placed at d = 1 Mpc, the UCHVCs are com- 
pact and do not violate any structure predictions of 
the ACDM scenario emphasized by SMW02. Their HI 
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HVC330.11 + • HI' i>'241 .0 + 1 1 5205 51 mJy km/s p be< 




Fig. 1.— Upper panel: Cloud at I = 330.11, b = +73.07, czq = 
124 km s^ 1 , W = 19 km s^ 1 , angular size of 6' X 4', integrated 
flux of 0.60 Jy km s" 1 , M HI = 1.4 X W b d{ Ipc Mq, M dyn [< 

Rhi] — 1.1 x 10 7 d,M pc Mq. Lower panel: HI line spectrum of the 
cloud shown above. The strong feature at zero velocity is the MW 
emission; the cloud is the narrow feature at czq = 124 km s _1 . 




Fig. 2. — Sky distribution of the UCHVCs (small diamonds, of 
area proportional to source flux. Galaxies with d < 2.6 Mpc are 
plotted as filled stars. A grid of Galactic coordinates is superim- 
posed. The star symbol outside the lower box is DDO 187. 

masses, < 10 6 Mq, would make them undetectable by 
previous HI surveys of nearby groups, such as that of 
Pisano et al. (2007). Thus the minihalo/HVC hypothe- 
sis could apply without posing the observational require- 
ment of LG uniqueness. HI surveys of nearby groups 
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Galactic Longitude 

Fig. 3. — Distribution of the UCHVCs in the Galactic longitude 
vs. Galactic Standard of rest plane, with superimposed lines of 
constant LSR velocities, 130 and -100 km s _1 , at two Galactic 
latitudes, b = 35° and b = 70°. Symbol code as in Figure 2. 

reaching a sensitivity level of 10 5 Mq — about an order 
of magnitude more sensitive than ALFALFA — would be 
necessary to verify the possible existence of similar sys- 
tems elsewhere than the LG. Such observations are cur- 
rently challenging but possible with foreseen upgrades to 
current instrumentation. 

The mean parameters of the UCHVCs at d = 1 Mpc 
are a good match for SMW02 minihalo models with a 
Burkert density profile, Rhi — 0.7 kpc, M HI ~3x 10 5 
Mq, total to neutral gas mass ratio of 15, peak Nhi — 
4 x 10 19 cm" 2 , total halo mass M mr ~ 3 x 10 8 Mq, 
surrounded by a hot, ionized IGM of presure Phim — 
10 cm -3 K,. The peak Nhi predicted by the models 
is higher than, but not incompatible with the observed 
values shown in Table 1, as the latter are averaged over 
the radius of the clouds and smeared by the ~ 3.5' beam 
of the Arecibo telescope. The Md yn inferred from the 
observations are about one order of magnitude smaller 
than the model's M V i r ; yet the two are consistent with 
each other since the values inferred from the observations 
are estimates of the masses within Rhi- As the cold 
baryons dissipatively collapse to the bottom of the halo 
potential well, the extent of the dark matter halo exceeds 
that of the cold gas by a factor of several, explaining the 
mass discrepancy. 

A further test of the minihalo/HVC hypothesis relates 
to the number of minihalo candidates expected from sim- 
ulations. The cumulative halo mass function N(> M), 
is a power law of approximate slope M~ 86 (e.g. Fig. 
6 of Gottloeber et al. 2003). Figure 4 of Hocft et al. 
(2006) shows the transition in the baryon fraction of ha- 
los dropping from the cosmic value of 0.16 to less than 
0.02 over a decade in halo mass, around the "characteris- 
tic mass" M c har — 10 9 ' 8 Mq, for which halos are able to 
retain half their baryons. If we assume that the UCHVCs 
pertain to the category of halos of mass between M c h ar 
and M c har /$m, and that the N opt optical galaxies in the 
field are hosted by halos with M > M c h ar , then we can 
estimate the ratio between the numbers of the two kinds 



<5^ 86 , with 5 



A I 



> 1. 



of objects as /„ = N opt /N c i oud 
The range in dynamical masses estimated for UCHVCs 
is about one order of magnitude (see Table 1); assuming 
that the corresponding, putative halo masses are spread 
just as narrowly, we guess 5m — 10. There are 5 galax- 
ies with d < 2.6 Mpc in the region in which the clouds 
were observed, and the exclusion of the spectral region 
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120 < Vq > —110 km s -1 in our search due to con- 
fusion with Galactic emission blanks about 50% of our 
search volume. The expected number of minihalos is then 
- 5 x 10 86 /2 - 18. We identified 27 clouds. Given the 
crudeness of the calculation, our cloud detection rate ap- 
pears compatible with a minihalo/HVC scenario. Next 
we ask: what else could the clouds be? 

In a Galactic fountain (Shapiro & Field 1976), gas is 
accelerated by strong stellar winds and supernova explo- 
sions in the disk and ejected to the Galactic halo to z- 
heights of a few kpc; gas clouds then cool and "rain back" 
onto the disk. At a distance of 3 kpc from the Galactic 
plane, the clouds' HI masses would be near solar, their 
sizes on order of 1 or few pc and crossing times less than 
1 Myr, much shorter than the ballistic timescales of 50- 
100 Myr. Our cloud sample is extracted fron the Galac- 
tic polar region, the component of the clouds' velocities 
perpendicular to the Galactic plane is dominant. Veloc- 
ities in excess of 200 km s _1 are difficult to accommo- 
date within a galactic fountain scenario with currently 
assumed values for the density of a galactic corona (e.g 
Fukugita & Peebles 2006; Grcevich & Putman 2009): 
small clouds would rapidly decelerate due to ram pres- 
sure and MW tidal forces would eventually disrupt them. 

Our survey region overlaps with the so-called "field 
of streams" (Belokurov et al. 2006), containing part of 
the tidally disrupted remnants of dwarf spheroidal MW 
satellites, namely the Sag and Orphan streams. With ve- 
locities between +20 and -100 km s _1 , their kinematics 
are very different from those of most of the UCHVCs, 
and an association is unlikely. 

The Magellanic Stream is witness of the tidal disrup- 
tion of the Magellanic Clouds by the MW. Most of the 
Stream is antipodal to the north Galactic polar region. 
However, its so-called "Leading Arm" (LA) extends to 
the northern Galactic hemisphere up to b ~ 35°, with 
velocities reaching +250 km s _1 . If the UCHVCs are 
a northern extension of the LA, according to the tidal 
model of Connors et al. (2006; see also McClure-Griffiths 
2008) they would be at distances of order of 100 kpc, have 
masses on order of 10 3 M Q and crossing times of <~ 10 
Myr. Tidal and ram pressure forces would be mild. This 
is a plausible scenario, although the LA would then ex- 
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tend twice as far and forward of the Magellanic Clouds 
than previous observations indicated. 

Galaxies grow via mergers and intergalactic gas infall. 
Infalling gas can be shock-heated to the virial temper- 
ature of the halo or fall at colder temperatures. Simu- 
lations show that the cold accretion mode tends to be 
more important for small mass halos: for a galaxy like 
the MW, > 90% of gas accretion is thought to take place 
in the hot mode (Keres et al. 2009; Dekel et al. 2009). 
Sancisi et al. (2008) discuss the evidence for extraplanar 
HI in several nearby galaxies and interpret it as resulting 
from a combination of mergers of gas rich satellites and 
cold gas accretion. In their HI maps, the extraplanar gas 
is seen at distances of 15 kpc or less from the main galaxy. 
Two difficulties arise with a model whereby our clouds 
are a manifestation of a similar phenomenon. First, the 
objection raised in the discussion of the galactic foun- 
tain scenario, that the clouds are too small, vulnerable 
and short-lived, holds in this case as well. Second, as 
shown by Mailer & Bullock (2004), thermal instability 
and conduction prevent the cooling of infalling gas on 
scales smaller than the so-called "Field length", which 
for the galactic corona translates to a mass of ~ 10 6 M . 
At perigalactic distances of ~ 100 kpc or less, our clouds 
would be much smaller than that limit. 

In conclusion, we report the discovery of a category of 
HVCs which are plausible minihalo candidates. They 
have properties and are found in numbers which are 
compatible with theoretical expectations for halos with 
masses < 10 9 M , and could not have been detected 
by extant HI surveys beyond the LG. However, it is not 
yet possible to exclude that they may be part of the 
wider scenario of the yet relatively poorly understood 
perigalactic HVC phenomenon. While difficulties arise 
in their interpretations with several more frequently in- 
voked models, reasonable adjustments could be made to 
fit some of them, most notably the possibility that they 
are an extension of the LA of the Magellanic Stream. We 
have thus discovered a category of objects consistent with 
the minihalo/HVC hypothesis, but we have not proved 
that such interpretation is unique. 
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